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ABSTRACT 

Context. Time-dependent gas-grain chemistry can help us understand the layered structure of species deposited onto the surface of 
grains during the lifetime of a protoplanetary disk. The history of trapping large quantities of carbon- and oxygen-bearing molecules 
onto the grains is especially significant for the formation of more complex (organic) molecules on the surface of grains. 
Aims. Among other processes, cosmic ray-induced UV photoprocesses can lead to the efficient formation of OH. Using a more 
accurate treatment of cosmic ray-gas interactions for disks, we obtain an increased cosmic ray-induced UV photon flux of 3.8 x 10 s 
photons cm^s" 1 for a cosmic-ray ionization rate of H2 value of 5xl0~ 17 s -1 (compared to previous estimates of 10 4 photons cm 4 s _l 
based on ISM dust properties). We explore the role of the enhanced OH abundance on the gas-grain chemistry in the midplane of 
the disk at 10 AU, which is a plausible location of comet formation. We focus on studying the formation/destruction pathways and 
timescales of the dominant chemical species. 

Methods. We solved the chemical rate equations based on a gas-grain chemical network and correcting for the enhanced cosmic 
ray-induced UV field. This field was estimated from an appropriate treatment of dust properties in a protoplanetary disk, as opposed 
to previous estimates that assume an ISM-like grain size distribution. We also explored the chemical effects of photodesorption of 
water ice into OH+H. 

Results. Near the end of the disk's lifetime our chemical model yields H2O, CO, CO2 and CH4 ice abundances at 10 AU (consistent 
with a midplane density of 10 10 cnr 3 and a temperature of 20 K) that are compatible with measurements of the chemical composition 
of cometary bodies for a [C/O] ratio of 0.16. This comparison puts constraints on the physical conditions in which comets were 
formed. 

Key words. Astrochemistry; Protoplanetary disks; Molecular processes; Comets: general; (ISM:) cosmic rays 



1. Introduction 

Understanding the evolution of gas-grain chemistry and the role 
of the dust size distribution in protoplanetary disks is paramount 
to understanding their history. This history is imprinted in 
the remnants of the formation of the disks, namely asteroids 
and cometary bodies in our solar system. However, in-depth 
theoretical and laboratory studies of gas-grain chemical net- 
work s have only bee n resea r ched in recent t imes dAikawa et aT 



I2009I) . which link the gas species to the icy surface of grains. 
We also use a refined approach to calc ulating the effect 
of cosmic rays on both gas and grains ( Leger et al] [l985; 



1997]: IWillacv et al] 1 19981: iRoberts et alJ 120071: IQbergetal 



2009; Visseretal. 2009; Woitk e et alJ 120091) . Chemical net- 



works in protoplanetary disk-like conditions have not been 
k> . as thoroughly studied as i n molecular cloud-like condition s 
j_| ; (lHasegawa & Herbstl 119931: iBatesi 119861: IWoodall et al.1 120071) . 
d ', For these reasons, we are interested in exploring chemical 
pathways in different regions of the disk, focusing on the role 
of OH, a highly reactive radical that can change the molecular 
composition of the gas when efficiently formed. To accomplish 
that, we study the chemical evolution and the chemical network 
in two regions of the disk midplane, including enhanced OH 
formation from cosmic rays and photodesorption of water ice. 

Our goal is to study the time-dependent chemistry in the 
cold comet-formation zone, which in a T Tauri-type proto- 
plane tary disk can be located at 10 AU dKellev & Wooden! 
120091) . We include in our chemical network the most important 
adsorption-desorption proc e sses for H2O, CO, CO2 and CH 4 
(lHasegawa & Herbstl 1 1993b ICollings et all 1200 4: Oberg etalJ 



ICecch i-Pestelli ni & Aielloll 19921: IShen et all2 004). Cosmic rays 
are important because they drive the chemistry in some regions 
of the disk via the ultraviolet (UV) field produced by their 
ionization of H^ dPrasad & Tarafdarll983blSternberg et alJl987t 
iGredel et al.l ll989). This secondary UV field can photodissoci- 
ate many chemical species and change the chemical balance, 
especially in regions of the disk that are opaque to stellar and 
interstellar radiation: the highly energetic cosmic rays can reach 
deep into the midplane of the disk. For this reason, the midplane 
of the disk can be considered as a cosmic -ray dominated region. 

To analyze the chemical evolution, the local environment 
of the chemical system should resemble the conditions in a 
protoplanetary disk as much as possible. One of the main 
tools for building chemical n etworks in astronom y, the Umist 
database for astrochemistry (IWoodall et al.l 120071) . lists some 
important chemical processes that have been calculated for 
molecular cloud-like conditions. The rates for some of these 
processes need to be recalculated for use in protoplanetary disk 
chemical networks, using the appropriate values for density, 
temperature, dust parameters, and radiation field. Cosmic ray- 
induced UV photo-dissociation is one example of the type of 
processes for which rates are calculated in the Umist literature 
using molecular cloud- l ike grain parameter s and H2 abundances 
dSternberg et alJ 119871: IGredel et a l. 1989). In protoplanetary 
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disks, dust grains are larger than in m olecular clouds due to 
aggregation and coagulation processes dD'Alessio et al . 2001; 
lNattaetal]|2007t iDominik etal] 120071) . Therefore, the overall 
grain surface and the local UV dust extinction will be smaller 
under protoplanetary disk conditions than under molecular 
cloud conditions. Here, we calculate the cosmic ray-induced 
UV (CRUV) photoprocess rates using the appropriate dust grain 
parameters. 



In the work by Dalg arno et al. (see ISternberg et al.l Il987t 
iGredel et al]fl987l fl989) on the interstellar medium and molec- 
ular clouds (on which the Umist06 dWoodall et alj|2007l) CRUV 
rates are based), the local extinction of cosmic ray-induced UV 
photons is dominated by dust grains; the gas contribution to 
the extinction is not considered to be important. However, in 
protoplanetary disks the gas composition can be very different 
from that of the interstellar medium or molecular clouds (e.g. 
high abundances of CH4) which can significantly increase the 
gas opacity. It is then possible that in some regions of the disk 
CRUV photons can be absorbed much more efficiently by the 
gas than by the dust. 

The general understand ing of photode sorption is getting 
more detailed n owadays jAndersson & van Dishoec kl 1200a 
lArasa et al]l20Toh . When a photon hits an adsorbed molecule 
it can photodissociate it, instead of directly desorbing the ice 
species into the gas phase. The products of this dissociation 
can either recombine on the surface or desorb individually. 
In the case of water, molecular dynamics studies of a pure 
water ice layer that is hit by a UV photon show that a water 
molecu le is often desorbed as OH+H i n stead of desorbin g 
intact ( Ander s son & van Dishoeckl 120081: lArasa et al.l l2010h . 
The effect of this constant OH formation channel through grains 
may drastically alter the gas-phase chemical evolution. 

We can safely neglect the gas opacity at distances of ap- 
proximately 10 AU in our model, because most of the material 
is frozen on the dust grain surface, thereby suppressing the gas 
phase abundances of many molecules. However, as we get closer 
to the star, the temperature is high enough to keep molecules in 
the gas phase. At 1 AU the gas opacity is therefore important 
for the chemical evolution. Another factor to consider is that 
gas-phase processes involving line absorption will be affected 
by the local temperature conditions. This implies that CRUV 
photodissociation cross sections will have to be recalculated 
and integrated properly into the rate equations. This will be the 
subject of a subsequent paper. 

The structure of this paper is the following. In Section [2] 
we explain how we use the r esults of the disk modeling code 
ProDiMo dWoitke et alj|2009l) as input for our model. Next we 
explain the role of cosmic rays (Section [3]), including an appro- 
priate calculation (i.e. following protoplanetary disk grain pa- 
rameters) of cosmic ray-induced UV photoprocesses in a low- 
gas environment. Section|4]consists of a comprehensive account 
of the gas-grain chemical model including adsorption-desorption 
mechanisms. Section|5]deals with the setup for the chemical evo- 
lution model, including an estimate of the initial set of chemical 
abundances that apply for a protoplanetary disk. In Section [6] 
we explore the full chemical networks at the comet formation 
zone, comparing the chemical abundances with those measured 
in cometary bodies. In Section|7]we discuss the chemical effects 
of H2O ice photodesorbing into OH+H instead of desorbing in- 
tact. A discussion of our results in Section [8] is followed by out 
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Fig. 1. Density structure (total hydrogen nuclei number density 
71(H)) of a T Tauri disk model as function of radial distance from 
the star r and the relative height z/r. Gas in thermal balance. 
The red dashed line shows the temperature boundary where 
T g - 1000 K. The regions under consideration for this paper are 
shown as black ovals (PI, and P2 for be nchmarking). This fi gure 
is taken from the ProDiMo simulation (Woi tke et al.l l2009). 



conclusions (Section|9]l, in which we summarize the main results 
of this study. 

2. Disk Model 

For the sake of consistency, the disk structure used here is 
taken from the model of a passive disk around a T Tauri star 
obtained with the steady-state disk modeling code ProDiMo 
dWoitke et al J 12009). This approach is valid because the accre- 
tion rate is negligible and transport processes occur on longer 
timescales than the ones modeled in this work. The position 
of the region under study within the disk structure (from the 
simulation) is illustrated in Fig. [TJ The location at 10 AU from 
the star corresponds to a likely location for comet formation, 
in which the density remains fairly high (« <H >=10 10 cm 4 ) if 
compared to molecular cloud conditions and the temperature is 
quite low (T-20 K). The CO and H2O adsorption-desorption 
processes are most significant at such low temperature^ 
Studying the change in chemical composition as a function of 
density while keeping a constant temperature is equivalent to 
moving in a direction perpendicular to the plane of the disk: disk 
models are vertically isothermal at high optical depth. Table [TJ 
summarizes all relevant input parameters. 

The grain size distribution in a quiescent protoplanetary 
disk is assumed to follow a power-law distrib ution f(a)~a~ p 
with p-3.5 and a in the range 0.1-10 /zm (Woitke et al. 2009). 
This value of p is also use d for grains in molecular clouds 
dWein gartner & Draindl200ll) . but Omin and a max are different, 
reflecting the overall smaller grain sizes in the interstellar 
medium (ISM): in the much denser protoplanetary disk en- 
vironment aggregation processes and dust settling to the 



In those regions, the gas a nd grain temperatures are coupled, i.e. 
< = Tdust ( Wo itke et all2 009). 
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midplane produ c e large r sized grains on average. For example, 
iD'Alessio et al.l (|2001) show the median Taurus SED and 
models it with an ISM grain size distribution and also with a 
distribution with larger grains on average. The latter is shown 
to fit much better. This grain growth causes a reduction of the 
average surface area of the grains, which limits their ability to 
absorb UV photons. 

Dust settling in the midplane could lea d to a lower gas-to - 
dust ratio than the typical value used here (Wo itke et alj 2009). 
This decreases dust extinction in the top layers in the disk 
(i.e. these layers become almost transparent), which allows the 
UV radiation to penetrate deeper towards the midplane. This 
means that the Ay -I line is shifted to a lower height above 
the midplane. However, this effect can be compensated by 
the increased dust UV extinction in the midplane due to the 
enhanced dust-to-mass ratio. To isolate the effect of larger dust 
grains on the CR-induced UV field, here we stick to the grain 
parameters used in previous ProDiMo papers and keep the value 
100 for the gas-to-dust mass ratio. 

The local UV flux (in erg-cm -3 ) of t he star and the ISM is 
frequently give n in units of a Draine field ( Weing artner & Draind 
l200ltlLee]ll984l) . The strength of the UV field is then provided 
via the dimensionless quantity ^ (which is 1 for the ISM), 
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The full 2D radiative transfer treatment in ProDiMo (which 
uses the same grain parameters, as stated in Table [TJ yields 
that for regions close to the midplane of the disk (z/r<0.05 and 
0.7<r<10 AU the disk is opaque to the stellar and interstellar 
UV photon s and thus the loca l UV field strength is very low, 
X - 0.001 (IWoitke et al.ll2009h . Hence, UV photoreactions will 
only play a minor role here when compared to the effects of 
cosmic rays. 

A cosmic ray-induced UV ( CRUV) field is present 
in most regions of the di sk dPrasad & Tarafdarl fl983t 
ICecchi-Pe stellini & Aielld [T992). This field can cause photodis- 
sociation reactions to take place at regions with high optical 
depth in the disk. However, CRUV photons are created locally 
and are assumed to be absorbed by the surrounding material 
(i.e. gas and dust grains) that is locally present. This means 
that the most appropriate treatment of CRUV photoreactions re- 
quires calculating the UV field created by cosmic-ray ionization 
of molecular hydrogen, accounting for dust grain surfaces as a 
possible sink for these UV photons by using the appropriate (lo- 
cal) grain properties, and also for the effects of gas opacity when 
the local conditions require. This will be discussed further in the 
following section. 



3. Cosmic rays 

Cosmic rays play a very important role in the chemistry be- 
cause they can pe rvade most regions of a protoplanetary disk. 
IShen et alJ (|2004) have calculated that a cosmic-ray iron nu- 
cleus with an energy ranging from 20-1000 MeV will lose less 
than 1 % of its energy when passing through a dense molecular 
cloud. Their influence can be intense in highly obscured regions 
(Ay > 5) where most stellar or interstellar UV radiation canno t 
penetrate dRoberts et al.ll2007l) . lUmebavashi & Nakanol d!981l) 
found that cosmic rays get effectively attenuated at a column 



Table 1. Table of modeling parameters, following IWoitke et al.l 
(120091) . The symbol # identifies an ice-phase species. 
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density above X ~ 150 g cm -2 . In the particular disk model used 
here, the column density at 1 AU is not high enough to provide 
a significant extinction of cosmic rays. Choosing a steeper sur- 
face density distribution or a smaller outer radius could make 
the midpla ne opaque to cosm ic rays at greater distances from 
the star. In Le ger et aD {l985) the total flux for cosmic-ray par- 
ticles with energies higher than 0.02 GeV/nucleon in the ISM is 
calculated to be 



^CR 



J0.02 



6 H (e) de ■ 



1.0 CR cm^ster's -1 



(2) 



The amount of cosmic raysthat are present can also be mea- 
su red from the CR ion ization rate of H2, £h 2 - A recent work 
by llndriolo et al.l (T2007) shows that column densities of H^ (a 
tracer of H2 ionization) vary in different molecular clouds. From 
their work, we pick a conservative value for £h 2 , which agrees 
with pr evious calculations obtained fro m measured cosmic ray 
spectra dCecchi-Pestellini & AieTloll 19921) 



f Hj =5x 10~" s" 



3.1. Impact on chemistry 



(3) 



Cosmic rays play an important role in molecula r cloud chemistry 
dHerbst & Klempererl [19731 IShen et al.l 120041) : reactions with 
the CR-generated H^ and He + are of special significance for 
CO/H2O gas formation and destruction in the gas phase. Both 
H^ and He + can dissociate CO and also SiO, which is impor- 
tant player in H2O formation. He + is created from He+CR reac- 
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tions, which happen at the same rate as H2 ionization (see Eq. 
(O) becau se they have appr oximately the same stopping power 
dHerbst & Klempererlll973l) . H+ is created after H 2 +CR colli- 
sions occur, which enable ionized molecular hydrogen to rapidly 
react with H2 



H+ + H 2 



H+ + H 



(4) 



As an important side effect, excited Hj will also emit Werner 
and Lyman UV photons that can photodissociate oth erwise 
stable species such as CH4 dPrasad & Tarafdarl 119831) . The 
presence of CR-induced UV photons can change the overall 
chemical balance through efficient ion-molecule chemical 
pathways. 

For processes such as photodesorption, we need to calculate 
the local CRUV flux, Nqu, while for CRUV photoprocesses we 
need to recalculate the efficiencies given the revised dust prop- 
erties. The effects of gas opacity can be neglected at 10 AU, as 
most of the material is frozen on the grain surface, and gas abun- 
dances are low. The next two sections discuss this in more detail. 



3.2. Cosmic ray-induced UV photon flux 

The effect of cosmic rays on molecular hydrogen is twofold: io n- 
ization and direct excitation (ICecchi-Pestellini & Aielld [l992) 



CR + H 2 -> e"(30eV) + H+ , 
CR + H 2 -> . 



(5) 
(6) 



Fifty-five percent of all electrons produced in reaction (|5]l 
will excite molecula r hydrogen to excited electronic states 
dSternberg et al.lll987l) : 



e +H 2 



H 2 . 



(7) 



Excited electronic states of H 2 via (|6]l and (Q) spontaneously de- 
cay to the B 'E* and the C 'E„ states, which then decay to the 
ground electronic state X 1 and emit Lyman and Werner UV 
photons in the wavelength range 90 - 170 nm: 



-» h 2 * + y uv . 

This is called the 
dPrasad & Tarafdarlll983l) . 



Prasad-Tarafdar 



(8) 

mechanism 



We calculate the flux of locally generated cosmic ray- 
induced UV photons (in cm" 2 s _1 ) based on the work of 
ICecchi-Pestellini & Aiellol 0992): 
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(9) 



Here u = Q S ca/(6sca + 2abs) is the grain albedo at 90 - 170 nm, 
and <x^ is the grain UV extinction cross section per hydrogen 
nucleus. The factor «h 2 /»<h> is included because CRUV photons 
are generated by H 2 molecules, so their number should depend 
on the local H 2 abundance. 

We calculate u> and cr|^ for the protoplanetary disk grain 
size distribution discussed in Section [2] using the appropriate 
dust opacity /cuv (see TableQ}: 



<H> 



Kuv- 



Pd 

«<h> 



1.47 x 10~ 22 cm 2 



(10) 



The Umist database uses a stan dard molecular cloud value for 
the UV extinction cross section ( Stern berg et alJlT987h 

- uv =2x 10- 



<H> 



-21 2 

cm 



(11) 



Such a divergence in <x|^ can be expected for a protoplane- 
tary disk-like grain size distribution, since it favors large grains 
that do not provide significant extinction at short wavelengths 
when compared to a molecula r cloud-like grain size distribution 
(Wei ngartner & Draind 12001). With this, it is possible to esti- 
mate a value for the CRUV photon flux under protoplanetary 
disk conditions, assuming that most of the H in the gas is stored 
in molecular hydrogeifl: 



Ncu - 386 000 photons cm""Y 



(12) 



Here Nm is 13 . 6 tim es more than the value implicitly used in 
ISternberg et al.l d 19871) (i.e. in the Umist database) and approxi- 
mately 40 times hig her than previous estimates for CRUV fields 
in molecular clouds (Shen et al. 2004): 



Nm 



10000 photons cirrV 1 



(13) 



Comparing our value for Nm with the stellar and interstel- 
lar UV field in the midplane (which is approximately 190000 
photons cm^s -1 ) and considering direct cosmic-ray ionization 
processes, it follows that cosmic-rays can be a major driver of 
the chemistry in the midplane of the disk. Aikawa et al. (1999) 
reached a simil ar conclusion based on the est imate of cosmic -ray 
extinction bv lUmebavashi & Nakanol d!981l) . 

3.3. Cosmic ray-induced UV photodissociation 

The secondary field of cosmic ray-induced UV photons dis- 
cussed in Section 13.21 has the same effect on the overall chem- 
istry as stellar or interstellar UV radiation, via photodissocia- 
tion/ionization of gas species and photodesorption of ice species. 
CR-induced photopr ocesses are liste d in the Umist dat abase 
dWoodall et alJ2007l) using data from lSternberg et alJ dl987l) and 
iGredel et al.l (Q~989). The cosmic-ray-induced photodissociation 
efficiency describes the competition for the locally generated UV 
photons between a species i and the rest of the gas and dust 
grains: 



„_ = r P(v)cr,-(v) 
J cr tot (v) 



dv. 



(14) 



Here cr,(v) is the photoprocess cross section dLedll984T) . P(v) 
the H 2 line emission probability (Lyman and Werner lines in a 
de-excitation transition of H 2 from a CR-induced 30 eV electron 
excitation), and cr tot (v) is the total extinction (of both gas and 
grain) cross section: 



o-tot(v) = o-Hy(l - oj) + V —<Tj(v) . 
w V «<H> 



(15) 



Dust UV extinction is assumed to be larg er than gas UV extinc- 
tion in molecular cloud-like con ditions ([Sternberg et al. 1987; 
ICecchi-Pestellini & Aiellol 1 1992b . However, for protoplanetary 
disk-type chemical abundances this only applies if the local gas- 
phase abundances are low, as is the case at 10 AU. Ignoring the 
contribution of the gas opacity, the efficiency is then written as 



f cr,(v)P,(v) 1 



(16) 



2 This remains valid as long as molecular cloud-like abundances are 
used as initial conditions of the protoplanetary disk model. 
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Table 2. Table of chemical species in the chemical networks. 



Type 


Symbol 


Atoms 


H, He, C, 0, S, Si, Mg, Fe 


Ions 


He + , Si + , Fe + , H~, H + , C + , 




+ , S + , Mg + 


Molecules 


H 2 , H 2 0, CH 2 , HCO, SiO, C0 2 , 




SiH, CH 3 , CH 4 , OH, 2 , CO, CH 


Molecular 


HCO+, CH+, H+, SiH + , SiO\ 


Ions 


CH+, H 3 + , H 3 + , SiHJ, CH+, 




CH+, H 2 + , SiOH + , CH + , H+, 




0+, C0 + , OH + , COJ 


Ice species" 


H 2 0#, CO#, CH 4 #, C0 2 # 



a: # stands for an ice-phase species. 



The reaction rate coefficient is usually calculated by considering 
that CRUV photons are generated in a proportional amount to 
the CR ionization rate of H2, 



1 



1 



(17) 



The previous expression ignores that CRUV photons are created 
from single H2 molecules, as opposed to hydrogen nuclei. For 
this reason the U mist06 database cor rects the photodissociation 
efficiencies from lGredel et al.l d!989l) by a factor 2, which is as- 
sumed to be the value for «(h>/»Ht in molecular clouds. For this 
reason we include the factor «Hi/«<h> explicitly in the CRUV 
photo-dissociation rate: 



kcv,i = 2 



& 2 



1 - u> \ 300 K 



»H 2 
»<H> 



(18) 



The exponent f3-\.\l accounts for the fact that some species 
are photo-issociated by discrete line absorption (IWoodall et alJ 
l2007t lGredel etal]ll987l) . For all other species other than CO, 
this parameter is zero. Using our new value for cr^ in Eq. 
(1101) . CRUV photodissociation processes will be up to 13.5 times 
more efficient due to the reduced grain extinctioqj. 



4. Chemical model 

The rate equations that describe the time-dependent chemistry 
(and chemical channels) of a gas-grain system comprise chemi- 
cal reactions of the following types: 

- gas-phase reactions: reactions between species that are 
present in the gas phase such as ion-molecule, recom- 
bination, charge transfer, neutral-neutral reactions, photo- 
dissociation and cosmic ray-induced processes. 

- adsorption-desorption processes: reactions that take species 
from the gas phase into the solid phase (on the surface of 
the grain) or vice versa, usually driven by the local temper- 
ature, cosmic rays, and the local UV radiation field (stellar, 
interstellar and cosmic ray-induced). 

- surface-surface reactions: reactions that occur on the surface 
of the grain, between species that are present in the solid 
phase. 



This factor is 13.5 when compared to the Umist rates, which are 
based on a standard molecular cloud-like dust UV extinction cross sec- 
tion. The factor 40 comes from the comparison to conservative esti- 
mates of the CRUV photon flux in molecular clouds, as in the previous 
section. 



Surface-surface reactions are beyond the current scope of this 
work, so only the first two are considered here. The 53 chemical 
species used here are listed in Table [2] 

The rate equation for the chemical network under study is a 
system of kinetically coupled differential equations. For a gas- 
phase species i, it is written as 

drij 
dt 

Here Rf is the rate that characterizes all processes that have the 
species i as a product, and Rj is the rate that describes all pro- 
cesses that decrease the abundance of the species i. For an ice 
speciefl /#, the rate equation accounts for adsorption-desorption 
processes: 



= RJ - RT 



(19) 



— — = IC; Hi - ki «,# 

at 



(20) 



The four most important desorption processes considered here 
are thermal, photo-, cosmic ray-induced photo-, and direct 
cosmic -ray desorption: 



kf = k m +k 



(21) 

It should be noted that we did not include X-ray desorption pro- 
cesses, which are very important near the inner rim, because the 
rates for these processes h ave not been determ ined beyond order- 
of-magnitude estimations (IWalshetalJl2010l) . 

4.1. Adsorption 

Adsorption is the process that allows a species in the gas phase 
to adhere to the surface of a grain upon collision. Since the grain 
is assumed to be immersed in a well-stirred gas at a temperature 
T„, a constant flux of gaseous molecules (each with mass mi) is 
hitting the surface of the grain with an average thermal velocity 
of 



(22) 



For a number density of dust grains «d us t, a sticking probability 
S — 1, and an average grain surface area of 47r(a 2 ) the rate is 



■ An{a 2 )S vfn dust 



(23) 



Adsorption becomes more efficient (i.e. acts on shorter 
timescales) as the gas density increases, given that more particles 
per volume are impinging on grains. The largest variation in the 
sticking paramet er S comes fr om the temperature-to-adsorptio n- 
energy ratio (see lFeglevll2000h . In lBurke & Hollenbachl(ll983l) a 
complete parameter-space exploration of the sticking probability 
yields that under protoplanetary disk-like conditions the chemi- 
cal balance does not change signif icantly, which see m to be con- 
firmed by experimental results in lOberg et al.l(l2009l) . Parameters 
for adsorption are given in Table [TJ 



4.2. Desorption 

Desorption processes counterbalance the complete depletion of 
species due to freeze-out. It should be noted that the details of 
the various desorption processes are not very well understood, 
are based on ve ry few laboratory e xperiments, and use extreme 
simplifications (Roberts et al. 2007|). 

4 The # symbol identifies an ice-phase (frozen on a grain surface) 
species /. 
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4.2.1. Thermal desorption 

Thermal desorption is a process in which the temperature of the 
dust grain can cause some of the ice-phase species bound to its 
surface to evaporate. The rate at which this happens depends on 
a characteristic binding energy 



fef th =v i exp 



~~kf~ dj 



(24) 



This Arrhenius-type expression for the thermal desorption rate is 
theoretical in nature, and the constant v, can be estimated from 
the vibrational frequency associated to the bond that holds the 
ice -phase species on the surface of the grain. This is based on 
the Polanyi-Wigner equation for a sing le desorption process, ig- 
noring rotational degree s of freedom dHollowav & Beebvlll975t 
lGalwev&Brownlll999h : 



V, = 



2n sad kEf 



Table 3. List of the most significant initial abundances of species 
used in the chemical networks. 



(25) 



Pre-molecular cloud 


Molecular cloud 


(Atomic) 


(Post-modeling) 


Symbol 


log(?ix/n(n}) 


Symbol 


log(nx/«<H>) 


H 





H 2 


-0.301 


He 


-1.125 


He 


-1.125 


O 


-3.538 


CO 


-3.939 


C 


-3.886 


H 2 0# 


-4.129 


Si 


-5.1 


O 


-4.56 


Fe 


-5.367 


co 2 # 


-4.575 


Mg 


-5.377 


o 2 


-4.977 


S 


-5.721 


SiO 


-5.1 






Fe 


-5.367 






Mg 


-5.377 






H 


-5.57 






H 2 


-5.798 " 






S 


-5.721 



Here n sur f = 1/A s i te = 1.5 x 10 15 cm~ 2 is the n umber densit y 
of available adsorption sites per unit grain area (Herbst 19931) . 
Parameters for thermal desorption are given in TableQ] 



4.2.2. Photodesorption 

Photodesorption occurs when a UV photon hits a specific ad- 
sorption site on the surface of the grain. The rate of this process 
can be written as 



7 d,ph / 2\ ^dust T7 t-t 

K i = 7X{a ) Ti^Draine . 

n act 



(26) 



This rate describes a Draine UV field^ impinging on the grain 
surface, and desorbing a species i. The photodesorption yield F, 
is measured in the laboratory u sing a UV field co mparable in 
strength to the Draine field (see Ober g et alj|2009l) . Given that 
the ice mantle on the grain surface is composed of many layers, 
the rate has to account for an incoming photon only being able to 
process a fraction of the icy molecules that are adsorbed onto the 
grain surface. This number of active sites on the grain surface is 



"act = 47r(a )n dust n surf 7Vi 



Lay 



(27) 



,d, P h 



4n SU rfN L 



-Yi. 



(28) 



Consequently, the rate for this process depends only on the 
amount of surface adsorption sites, the photodesorption yield, 
and the local flux of UV photons. This process can become 
important when the timescale for thermal desorption is large, 
such as for H2O at temperatures below 120 K. Parameters for 
photodesorption are given in Table Q] 

Most species are assumed to photodesorb intact. However, 
this is not always the case for water, as reported i n mo lecu- 
lar dynamics studies by lAndersson & van Dishoeckl d2008h and 



5 A valu e for the flux produc ed by this field Forainc - 2 x 10 s cirr 2 s~' 
is given in Woit ke et alj J20 09). The wavelength range for this field is 
given in Section|2]above. 



lArasa etal ] (120101) . Most of the time (70%) a water ice molecule 
that is hit by a UV photon will dissociate and subsequently des- 
orb into an OH molecule and an H atom. This means that there 
is only a 30% chance that the water molecule will desorb in- 
tact. It should be noted that this special case of water ice desorp- 
tion does not change the photodesorption yield value. Here it is 
assumed that desorption occurs only for the first two layers of 
frozen material on the grain surface. When considering deeper 
layers, less than half of the water ice molecules desorb as water 
vapor. The effect on the overall chemical evolution of this special 
desorption case for water will be explored in Section|7] 

4.2.3. Cosmic ray-induced U V photodesorption 

This photodesorption process is identical to regular photodes- 
orption, except for the source of UV photons. Therefore, by re- 
placing the stellar and interstellar UV energy density ^^Draine in 
Eq. ( 1281 ) by the appropriate CRUV photon flux for a protoplane- 
tary disk Ncu in Eq. i 121) . the cosmic ray-induced UV photodes- 
orption rate becomes 



Here N^y = 2 is the number of layers of ice that can be affected £ d 
by an incoming UV photon. The expression for the photodesorp- ' 
tion rate in Eqn. ( l26l i can then be rewritten as 



4n su rfA ? L E 



(29) 



ay 



4.2.4. Direct cosmic ray desorption 

Grain heating by 



cosmic 



rays 



described by 

lHasegawa & Herbsil JT993) happens when a cosmic ray 
passes through a disk and hits a dust grain, heating it up to 
an estimated temperature of 70 K. This causes species to 
thermally desorb from the grain as they would from thermal 
desoption at 70 K. Iron nuclei have bee n found to be the most 
important contributors to grain heating (lLegeretal] H985). For 
an iron-to-hydrogen ratio of approximately [Fe/H]^1.6xl0~ 4 , 
it is possible to estimate the rate at which cosmic rays hit an 
a ^ 0.1 /jm grain: 



R CR = d>£ e R 7r<a 2 ) = 3.16 x 10" 14 s _1 



(30) 



The timescale for succesive CR-grain hits is then tqr — 10 6 yr. 
Given that the estimated cooling time for a silicate grain that 
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goes from T-7Q K to T-2Q K is on the order of 10 5 s, we can 
find the ratio between the cooling and heating timescale^] with 

/(70K) = — = 3.16 x 10~ 19 . (31) 

TCR 

This is called the "duty-cycle" o f the grain heating by cosmic 
rays (Hasegawa & Herbsn IT993I) . The "duty-cycle" modulates 
the thermal desorption at T-70 K, scaled to the CR ionization 
rate of H2 

£?' cr = /(70K)flf th (70KK H2 . (32) 

5. Gas-grain chemistry model 

The time-dependent chemical network used for this model 
is chem_compact, which is base d on Milica Milosa vlevic 
& Inga Kamp's chemistry code dMilosavlevid 1200 8). This 
code is a time-dependent solver of the gas-phase chemical 
rate equation and has been originally used to explore the role 
of shocks in protop lanetary disks. It use s the Umist database 
for astrochemistry (IWoodall et alJ 120071) . which identifies all 
relevant chemical reactions for a set of chosen species. We use 
the ordinary differential equation (ODE) solver vode, which is 
specifically designed to solve stiff ODEs with strong and sudden 
time variati ons, such as the rat e equation for a chemical network 
of any size dBrown et alJI 1989T> . This code has previously been 
benchmarked against steady state abundances, an d in addition 
we m ake our own benchmark against ProDiMo (Woitk e et alJ 
2009) steady state abundances in Section l572l 

Additionally, we include adsorption-desorption reactions for 
CO, CO2, CH4, and H2O, and a new calculation of grain pa- 
rameters for an appropriate treatment of CRUV photoprocesses 
for protoplanetary disk-like grain parameters (see Section [3]). 
We checked for the effects of incorporating other significant ice 
species such as O2 and SiO ice to our network, and found no sig- 
nificant differences from our results. At 10 AU, O2 and SiO ice 
are formed, but because adsorption is very efficient, it is neces- 
sary for them to be formed in the gas phase first. In future work, 
when we also consider surface reactions, we will incorporate ad- 
sorption for all gas-phase species. 

5.1. Initial conditions 

Assuming that the material has been processed before the for- 
mation of the disk, the initial abundances for the protoplanetary 
disk model were computed from molecular cloud-like condi- 
tionfl The chemical model was evolved from atomic conditions 
(see Table [3] left column) in a gas of density n<H> = 10 6 cm -3 , 
temperature T = 20 K and x = 0.01 for 10 7 years. The final 
chemical abundances for this run (see Table [3] right column) 
were used as input for the initial chemical abundance^] in the 
different protoplanetary disk models studied here. A similar 

6 Calculation of this value requires knowledge of the specific heat 
obtained from indirect measurements, which means that it is highly de- 
pendent on T and thus may not be valid for high (> 150 K) dust tem- 
peratures. 

7 This means that the molecular cloud values for the UV grain extinc- 
tion cross section (Eq. lilt were used in the CRUV photo-rates (Eqns. 
[Tj|]and|29}. 

8 The O2 and H 2 abundances under molecular clou d-like conditions 
are in consistent with current observations of the ISM (Hollenbach et al. 
2009). However, we found that formation of O2 and H 2 in our proto- 
planetary disk model does not depend on their initial abundance. 



approach has been used bv lThi et"aTI (1201 ll) . 

We studied the chemical evolution under two different con- 
ditions for each point of interest: using a low (molecular cloud- 
like) and high (protoplanetary disk-like) value for the CRUV 
field. Thus we can see how the chemical abundances and path- 
ways change as the protoplanetary disk-sized grains absorb less 
CRUV photons. 

5.2. Chemistry benchmarking at 1 AU 

For the sake of consistency, we teste d the solver again st the 
steady-state chemistry in ProDiMo (IWoifke et all I2009h . by 
comparing the results at a distance of 1 AU from the star, near 
the midplane of the disk (P2 in Fig.[Q: T = 80 K, n <H > = 10 14 
cm -3 and^ = 0.001. For the purpose of matching the results of 
both models, we used a low (molecular cloud-like) value for the 
cr^ and ignore the effects of CRUV desorption. 

In the PRoDiMo-simulated chemistry at 1 AU all the oxygen 
is trapped in FbO ice («h,o#/ w <h> - 10~ 4 ), and all the carbon in 
methane gas (ncH 4 /"<H> — 10~ 4 ), which implies a low abundance 
of gas-phase CO («co/«<h> — 10~ 6 ) and consequently a very low 
atomic oxygen abundance («o/«<h> - 10~ 12 ). chem_compact 
yields very similar results, although on a much larger timescale 
than the lifetime of the disk (r ~ 10 8 yr). Therefore in this case, 
steady state abundances in the midplan e have to be t reated with 
caution. This has been pointed out by I Woitke et alJ |2009)- see 
their Fig. 13. 



6. Chemistry in the comet-formation zone 

In the midplane of the disk at a distance of 10 AU from the 
central star, the density of our model is m<h> = 10 10 crrT 3 and 
the temperature is T — 20 K, corresponding to the region PI 
in Fig. Q] Even though photo-desorption here is more efficient 
than thermal desorption, it acts at a timescale much longer than 
the disk's lifetime. These conditions are ideal for ice formation 
on grains, which means that carbon and oxygen will be trapped 
on the surface of grains and will not efficiently form gas-phase 
molecules. The chemical evolution in this region of the disk 
is plotted in Figs. [2] and [3] corresponding to a low and a high 
CRUV flux, respectively. 

From these figures it follows that while the local CRUV 
field does not radically change the chemical balance, it does 
have an effect. This will be discussed in Section 16721 It should 
also be noted that cosmic -ray ionization of H2 and He has a very 
big impact on the chemistry, even though those rates are not 
affected by the dust grain parameters. The formation pathways 
that we describe in the following section apply to both low and 
high CRUV flux conditions. 



6.1. Chemical pathways 

The final abundances and time evolution of CO shows that it 
suffers depletion on a typical timescale of 10 6 yr. This suggests 
that H2O is more efficiently formed than CO on long timescales, 
removing oxygen from the gas phase and trapping it into water 
ice. The main processes that create CO, H2O and other related 
species are illustrated in Figs.|4]and|5] where the gray pathways 
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1 2 3 4 5 6 7 

log(time/yr) 



Fig. 2. Example of the time evolution of gas/ice abundances 
(thin/thick lines) using the standard ISM value for cr^ — 2 x 

10~ 21 cm 2 . The disk parameters correspond to point PI in Fig.Q] 
(n <H > = 10 10 cm" 3 and T = 20 K). 
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Fig. 3. Example of the time evolution of gas/ice abundances 
(thin/thick lines) using an appropriate <x^ = 1.5 x 1CT 22 cm 2 
for a protoplanetary disk. See text in Section ROl for a discussion 
of these results. The disk parameters correspond to point PI in 
Fig.Q](n<H> = 10 10 cirT 3 and T = 20 K). 



indicate CRUV photodissociation. 

Figure [3] shows that after 100 yr, CO ice is the main carbon 
carrier. Excluding adsorption, the main reaction destroying gas- 
phase CO is 



CO + He 4 



C + + O + He . 



(33) 



It should be remembered that He + is created by cosmic -ray ion- 
ization of He. C + created in this reaction will undergo charge 
exchange with other atoms (Si, Mg, Fe) to form atomic carbon, 
which then reacts very efficiently via radiative association with 
H 2 toformCH 2 : 

C + H 2 -> CH 2 + y ■ (34) 

In turn, CH 2 reacts with atomic oxygen to again form CO: 

CH 2 + O -> CO + H 2 . (35) 

This reaction is very fast because atomic oxygen is steadily pro- 
duced by He + dissociation of CO in reaction ( 133b . CO is also 
formed from HCO + , either in the Si-reaction, 



HCO + + Si -> CO + SiH + , 
or via HCO 

HCO + + Fe, Mg HCO + Fe + , Mg + , 
HCO + O -> CO + OH , 
an HCO + can itself be formed from CO via H 4 ; 
CO + H 4 -> HCO 4 + H 2 . 



(36) 

(37) 
(38) 

(39) 



However, this HCO 4 -CO feedback cycle is not closed because 
C0 2 can also be formed from HCO: 



HCO + O -> C0 2 + H 



(40) 



After formation, C0 2 is rapidly frozen, at the expense of a 
significant fraction of CO. This reduces the carbon available for 
CO ice formation. 



C0 2 ice is more abundant in Fig. [3] than in Fig. [2] at the ex- 
pense of water, HCO and CH4. The long timescale for CH4 ice 
formation is due to the chemical pathways that create CH4 from 
CO. They are activated via CR-induced dissociation of CH3, 
which acts on a timescale of approximately 8 x 10 6 yr. CH3 is 
created from HCO and CH 2 



CH 2 + HCO -> CO + CH 3 



(41) 



The first and most important pathway starts from H 4 /- 
dissociation of CH3. It should be noted that H 4 / is created after 
cosmic-ray ionization of H 2 : 



CH 3 + H3 — > H 2 + CH 4 (+H 2 -> CH 4 + H) 



(42) 



The second pathway is started by CRUV photo-ionization of 
CH3 (which is half as efficient as He 4 dissociation): 



CH 3 + rcu -> e" + CH 4 , (+H 2 -» CH 5 4 ) 



(43) 



Thus, if CH3 is efficiently dissociated, CH 4 will be efficiently 
formed. CHt is only one step away from forming CH4: 



CH 4 + CO -> HCO 4 + CH 4 



(44) 



This reaction is more efficient at forming HCO 4 than reaction 
(l39l . After CH4 is formed, it is rapidly adsorbed onto the grain 
surface. 

SiO forms H3O 4 in reactions with He 4 , which can form 
SiO again via OH. This feedback cycle is broken on the same 
timescale as the freeze-out of water, which causes the SiO abun- 
dance to decrease on long timescales (see Figs. [2] and [3]). The 
SiO/H 2 feedback cycle goes as follows: 



SiO + He 4 -> 4 + Si + He , 
O 4 + H 2 -> OH 4 + H , 
OH 4 + H 2 -> H 2 4 + H , 
H 2 + + H 2 -> H3O 4 + H . 



(45) 
(46) 
(47) 
(48) 
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Fig. 4. Final (10 7 yr) chemical network for CO and CH 4 at PI 
in Fig. Q] (n <H > = 10 ll) cm" 3 and T - 20 K), where the low 
temperature freezes CH4 after its gas-phase formation at 10 6 
yr, as seen in Fig. [3] The lines represent reactions between 
species on the diagram, with arrows pointing at the product(s) 
of each reaction. Gray lines: Cosmic ray-induced photodissocia- 



tion reactions. (Notes: a. Via ycRuv/Hj 
HCO + +O^C0 2 +H or HCO + +O^CO+OH.) 



CH+/CH 



CH+ b. 




H 3 + 


He+ 


SiO 







0+, OH+, H,0+ 





Fig. 5. Final (10 7 yr) chemical network for H2O at PI in Fig.Q] 
(«(H> = 10 10 cm 3 and T = 20 K). The lines represent reac- 
tions between species on the diagram, with arrows pointing at 
the product(s) of each reaction. Gray lines: Cosmic ray-induced 
photo-dissociation reactions. 



Then, H 3 Q + efficiently forms H2O via dissociative recombina- 
tions dBateslll986l:ISternberg et al.lll987l) : 

H 3 + + e — > H + H z O . (49) 

This reaction creates water vapor fairly efficiently on a timescale 
of 0.06 yr. However, there is another dissociative recombination 
reaction with Si, driven by the high abundance of atomic silicon 
formed in reaction (f43T >: 

H 3 + + Si -» SiH + + H z O . (50) 

The previous cycle continously generates water vapor that 
can rapidly be adsorbed onto the grain surface. This chain 
of reactions is very efficient because no other processes are 
creating + , OH + or H 2 + . 

On the other hand, water (in the gas phase) can go back to 
H 3 + via H+ 

H 2 + Hj -> H 3 + + H 2 . (51) 

The dissociative recombination of H 3 + and an electron can also 
create OH 

H 3 + + e - -> OH + H 2 . (52) 

OH is a low-abundance but rapid catalyst for SiO formation 

OH + Si -> SiO + H . (53) 

However, since OH forms SiO about ten times faster than SiO 
forms H 3 + , SiO will be steadily formed within 10 6 yr, despite 
constant dissociation by CRUV photons. When we freeze all 
the water formed in reaction d50l >. then the reaction ( Bil l will 
be interrupted, and the OH-fueled SiO formation cycle is bro- 
ken. Thus, even though CRUV photodissociation of water va- 
por into OH and H can favor SiO formation over water on long 
timescales, in neither case (low or high CRUV field) can we 



produce enough OH to keep high abundances of SiO in the gas 
phase after a few Myr, as seen in Figs.[2]and[3] SiO is depleted 
after a few Myr because its formation pathway, reactions d45l l to 
(l50l . is disrupted by freeze-out of water. 

6.2. The effect of CRUV enhancement 

Besides being important drivers for the later evolution of the 
chemistry, as seen in the previous section, CRUV photons 
affect different species in different ways. This is evident when 
comparing the chemistry that arises in the low and high CRUV 
flux environments, as seen in Figs. [2] and [3] For instance, since 
CH 3 has an estimated higher cross section than CO and CH4, it 
is more susceptible to CRUV photodissociation. Thus, when the 
CRUV flux is enhanced, the CH4 formation pathways will be 
slowed down and its abundance will decrease. Also, even though 
there will more carbon available for CO formation, gas-phase 
CO will be more efficiently dissociated than C0 2 . This shows in 
Fig. [3] where we see the (final) C0 2 ice abundance increasing 
more (27%) than the gas-phase CO (7%) in the low CRUV flux 
case, in Fig. [2] 

This C0 2 ice enhancement comes at the expense of water, 
which decreases about 20% in abundance because of CRUV 
photodissociation . This process enhances the OH abundance, 
and this shows how SiO can survive for a few more million years 
(as noted in the previous section). This shift in abundances is 
further proven by the long-term availability of atomic oxygen 
(formed in the photodissociation of CO and SiO), which in the 
high CRUV flux case is depleted a few hundred thousand years 
later than in the low CRUV flux case. 

7. Alternative H 2 photodesorption mechanism 

Molecular d ynamics sim ulations by Anders son & van Dishoeckl 
(12008b and lArasaetafl d2010l) show that the effect of UV 
photons impinging on water ice molecules is more destructive 
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Fig. 7. Final (10 7 yr) abundances of significant species at the 
comet-forming region (PI in Fig.Q]). Triangles: Water ice pho- 
todesorbs intact (Fig. [3}. Squares: Water vapor is partially pho- 
todissociated into OH + H. Plot obtained using an appropriate 
cr^ for a protoplanetary disk and for a [C/O] ratio of 0.46. 



Fig. 6. Schematic chemical network for a protoplanetary disk, 
showing the competition of some chemical species (H2O, CO, 
CH3 4, SiO, O2) for the carbon, silicon, and oxygen, depending 
on whether there is high OH formation or not. The lines rep- 
resent reactions between species on the diagram, with arrows 
pointing at the product(s) of each reaction. The dotted branches 
correspond to a chemical network where OH is not being ef- 
ficiently formed. The dashed branches replace the dotted path- 
ways when OH is being efficiently created, thus taking away 
oxygen from H2O to form CO, SiO, O2, and atomic oxygen. 
The gray arrows represent CRUV photodissociation processes. 



than previously thought. Before these studies, water molecules 
were believed to desorb intact into the gas phase after being 
expelled from the grain surface by a UV photon. However, 
water molecules are not desorbed instantaneously after the UV 
photon hits them: almost every single photon-water molecule 
interaction will result in photodissociation of the water molecule 
into OH and H. Most of the time (70%), these byproducts will 
leave the surface in the gas phase. However, it is also possible 
(30%) that the OH and H recombine on the surface, and the 
energy left over from this reaction will cause the newly formed 
water ice molecule to desorb into the gas phase. 

The main processes that create CO, H2O, and other related 
species are summarized in Fig. [6] schematically showing the 
chemical network depicted in Figs. [4] (for CO, CO2, and CH4) 
and [5] (for H2O and SiO). They also show CRUV photodisso- 
ciation processes and the new pathways arising from having an 
efficient OH formation mechanism. 

When water ice photodesorbs intact at PI (Fig. [7] corresponding 
to the final abundances in Fig.0 most of the oxygen is in H2O, 
CO, and CO2 ice, and a fraction of carbon is in CH4 ice (see 
O. When water ice partially desorbs into OH and H (Fig. [7]l 
the efficient OH formation enables the oxygen in water to be 
stored elsewhere, such as in atomic oxygen and SiO. Now all 
the carbon is trapped in CO and CO2 thanks to OH-driven CO 



formation (see Fig. [6]). 

The main chemical products of a high CRUV field present 
are atomic oxygen and, to a lesser degree, OH. Atomic oxygen, 
which is now very abundant will react very efficiently with des- 
orbed OH to form O2 



O + OH -> 2 + H 



(54) 



This means that carbon and silicon photodissociated from CO 
and SiO will react with O2 to form CO and SiO at approximately 
the same rate (while other carbon carriers like CH4 are rapidly 
photo-dissociated): 



c + o 2 -» CO + o , 



Si + Q 2 



SiO +0 



(55) 



(56) 



This feedback pathway arising from a highly efficient OH- 
forming mechanism enables the long-term gas-phase formation 
of O2 and SiO that can later freeze on the dust grain surface. Now 
that carbon forms CO via O2 and not via CH2 as in reaction ( l35l >. 
the byproducts CH3 and CH4 are not being efficiently formed, 
and CO becomes the main carbon carrier. The final chemical 
abundances resulting from having this special case of water des- 
orption in PI (see Fig.[TJ is depicted in Fig. [7] 



8. Discussion 

In this work we focus on the timescales in which gas-grain chem- 
ical processes can take place in passive, nonaccreting disks. We 
thus assume that the density and temperature of the midplane 
will stay approximately constant during our simulation. We dis- 
cuss the main results drawn from our study of the gas-grain 
chemical timescales in relation to observables such as the com- 
position of cometary ices. 



10 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



8.1. Implications for comet formation 

Ice formation for different molecules takes place at different 
times because desorption processes act different timescales. This 
is evident in Fig. [3] where CO ice is more abundant than water 
ice between 10 2 and 10 6 yrs. After this t he ice content resembles 
what is observed in cometar y ices0 (Bockelee-M orvan et alj 
2004; iBockelee-Mor van 2010), where the predominant form of 
ice is H2O, followed by CO and CO2 ice and, to a lesser degree, 
CH4 ices. CH4 ice is formed after a few million years, and only 
after water ice formation has become so efficient that it takes 
away most of the oxygen from CO and leaves the carbon free to 
form other molecules. This suggests that the CH4 ice found in 
comets was formed very late in the disk evolution. 

A comparison between measured cometary ice abundances 
and our resulting ice abundances is presented in Table [4] Our 
resulting «co#/"co,# and «ch 4 #/«h 2 o# ratios are very close to 
the cometary values, which means that the relative abundances 
of ice species are predicted by our model to some extent. 
Furthermore, when we attempted to increase T (i.e. moving 
radially toward the star) or reduce w<h> (i-e- moving vertically 
away from the midplane) in order to decrease the adsorption 
rates, we ended up with a very different n<x>#/ n co 2 # ratio, due 
to the high volatility of CO compared to CO2. This means that 
the conditions at PI are optimal for formation of comet-like ice 
abundances. 

This conclusion is unique to the density-temperature com- 
bination rather than to the exact location in the midplane of 
the disk. If the temperature changes by 10 K, different ice 
abundances will vary in a nonuniform way due to the different 
binding energies of the chemical species; for example, at 10 
K CH4 ice is about 100 times less abundant because it cannot 
be formed efficiently. At 30 K thermal desorption of CO is 
so efficient that most of it stays in the gas phase. Also, if the 
density changes, the gas chemistry that drives the formation 
of ices via adsorption changes, so the particular ice structure 
will also change. At a density of 10 11 cm 3 , CO adsorption is 
so efficient that the gas-phase CH4 formation processes are not 
efficient, and again lead to a CH4 ice abundance that is 100 
times less than in the 10 AU case. At a density of 10 9 cirT 3 , CO2 
ice becomes twice as abundant as CO ice. Thus, a difference of 
an order of magnitude in density or 10 K in temperature will 
drastically change the ice composition. 

The main reason for the discrepancy betw een the measure- 
ments compiled by Boc kelee-Morva 3 J20T0I) and our results 
stems from the fact that the carbon-to-oxygen ratio is much 
lower in observed cometary ices (~0.16) than in our model 
(~0.45). The reason behind this could be that some carbon is 
bound in dust that we consider here to be in the gas phase. In 
other words, as the sum of the carbon in the ices and in the gas 
phase is the total amount of carbon considered in our model, a 
fraction of it may be stored in refractory grain cores. Because 
of this we changed the [C/O] ratio in our model to the cometary 
value, which yielded the results in the right column of Table |4] 
These results fit the cometary nco#/«H,o# and «co,#/«h 2 o# ratios 
more closely. 



Table 4. Ice ratios in comets compared to our results at two dif- 
ferent epochs. 



Ratio 


Cometary 
measurements" 


c.c*: [C/O] =0.45 
10 6 yr 10 7 yr 


c.C : [C/O]=0.16 
10 7 yr 


n C0 2 #/ M CO# 


<0.6 


0.65 1 


0.43 


«CO#/«H 2 0# 


<0.2 


0.6 0.45 


0.14 


WC0 2 #/«H 2 0# 


0.02-0.12 


0.39 0.6 


0.06 


ncH 4 #/«H 2 o# 


0.003-0.015 


0.001 0.06 


<0.001 



a: Abundances compiled f rom cometary measur ements in 
IBockelee-Morvan et alj d2004h : iBockelee-Morvanl d2010h . b: Ratios 
from our chem_compact (c_c) model using a [C/O] ratio of 0.45 (see 
Fig. [3]l. a: Ratios from our chem_compact model using a (cometary) 
[C/O] ratio of 0.16. 



9. Conclusions 

Our time-dependent chemical networ k is compiled from 
the Umist dabase for astrochemistry dWoodall et al.l 120071) . 
adsorption-desorptio n processes for H2O, CO, CO2, and CH4 
dWoitke et alj l2009h. and an appropriate treatmen t of CRUV 
photoprocesses ( Cecchi-Pestellini & Aiellol 1 19921) for proto- 
planetary disks. 



9 Even though ice ratios can vary among individual comets (see 
Bockelee-Morvan et al. 2004), the general composition «h 2 o# > w co# > 
«co 2 # > «ch 4 # is fairly consistent. 



In the midplane of a protoplanetary disks, the steady for- 
mation of H2O and CO and their relation to secondary oxygen 
and carbon carriers, such as CO2, CH3, CH4 and SiO in the 
gas phase, are caused by recombination of ionized species. The 
presence of ionized material comes from an internal UV field, 
which in the dark, cold midplane reaches the floor level caused 
by cosmic ray interactions with the gas. Thus CRUV photons 
become the main driver of the chemistry. Locally generated UV 
photons can either photodissociate a molecule in the gas phase 
or impinge on the surface of a dust grain and photodesorb a 
frozen molecule. 

By calculating grain parameters (such as the UV albedo and 
extinction cross section) for the grain size distribution appropri- 
ate for a protoplanetary disk, we find a CRUV photon flux of 
380000 photon cm -2 s -1 , which is 40 times larger than conser- 
vative estimates for molecular c l ouds ( Prasad & Tarafdar 1983; 
I Cecchi-Pestellini & Aiellolll992l;ls"hen et al.ll2004l;lRoberts et al.l 
2007). Also, CRUV photodissociation processes can be up to 
13.5 times more efficient in the midplane of a protoplanetary 
disk than in a similar molecular cloud-like environment. This 
can be curbed by the effects of gas opacity in warmer regions 
of the disk (closer to the star than 10 AU), which will be the 
subject of our next paper. 

Cosmic ray-induced UV photons are responsible for the 
destruction of CH3 (which favors CO formation) and formation 
of CH4. Since CH4 and CO2 are more sensitive to CRUV 
photodissociation than CO, CO will mostly compete with H2O 
as the most abundant oxygen-bearing species for the region of 
the disk that we probed here (r ~ 10 AU). Carbon and oxygen 
are efficiently trapped in ices in the midplane of the disk because 
the desorption timescales are longer than the lifetime of the disk. 

The chemical abundances and their evolution change drasti- 
cally whenever there is a mechanism that efficiently forms OH. 
A high formation rate of OH implies that O2, SiO, and atomic 
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oxygen will be efficiently formed, often at the expense of part 
of the oxygen in H2O. Also, all the carbon will tend to be stored 
in CO instead of CH4. 

At 10 AU the only way to enhance OH formation is by pho- 
todesorbing water into OH+H. Even an enhanced CRUV field 
does not have any effect on OH formation because adsorption 
of gas molecules as soon as they are formed is very efficient at 
20 K. When OH is not formed via photo-processing of water 
ice, we obtain ice-on-grain abundances and ratios that are com- 
parab l e to those measured in comets ( Bockele e-Morvan et alJ 
2004; Bockele e-Morvanl |20 1 Ol) . One important instance is the 
formation of CH4 ice, for which we obtained a timescale of a 
few Myr. This timescale is limited by (cosmic-ray generated) 
He + dissociation of CH3. 

The sensitivity of these ice ratios to temperature and density 
provide strong evidence that cometary ices must have formed 
under conditions similar to the ones used here, and most likely 
at late evolutionary stages (after a few Myr). This also means 
that the measured ice composition of comets precludes their for- 
mation in an OH-rich environment. 
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